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Abstract. The Keller-Segel model is a system of partial differential equations 
modelling chcmotactic aggregation in cellular systems. This model has blowing up 
solutions for large enough initial conditions in dimensions d > 2, but all the solutions 
are regular in one dimension; a mathematical fact that crucially affects the patterns 
that can form in the biological system. One of the strongest assumptions of the Keller- 
Segel model is the diffusive character of the cellular motion, known to be false in many 
situations. We extend this model to such situations in which the cellular dispersal 
is better modelled by a fractional operator. We analyze this fractional Keller-Segel 
model and find that all solutions are again globally bounded in time in one dimension. 
This fact shows the robustness of the main biological conclusions obtained from the 
Keller-Segel model. 
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1. Introduction 

Biological pattern formation is a topic of growing interest in the field of applied 
mathematics, both due to the possibility of developing new mathematics as well as 
for the broad range of important applications it might have pLj. At the cellular level, 
it is known that chemotaxis plays a fundamental role in the self-organization of many 
biological systems. Chemotaxis is the directed movement of an organism in response 
to ambient chemical gradients, that are oftenly segregated by the cells themselves. In 
those cases where the chemical products are attractive (and they are therefore called 
chemoattractants), they lead to the phenomenon known as chemotactic aggregation: 
the cells accumulate in small regions of space giving rise to high density configurations. 
The patterns appear therefore as strong density variations in the spatial distribution of 
the cells. 

One of the most important partial differential systems for understanding 
chemotactic aggregation is Keller-Segel model [2], which in a simplified form reads: 



Here Dp is the cellular diffusion constant, k the chemotactic coefficient, 7 the rate 
of attractant production, (3 the rate of attractant depletion, Dc the chemical diffusion 
constant, p is the cell density, and c is the chemical density. The terms in Eq. fllap include 
the diffusion of the cells and chemotactic drift. Eq.(fT6l) expresses the diffusion and 
production of attractant. This system is known to have finite time blowing up solutions 
for large enough initial conditions in dimensions d > 2, but all the solutions are regular 
for d = 1; there is a large literature on the analysis of this problem for the Keller-Segel 
model as well as for some simplifications [31 HI [5l El [TJ [H [U [TOl [TTl [12] . The biological 
meaning of this mathematical fact was examined in Refs. [TH [H] . One conclussion of 
these works is that in a three-dimensional system, while collapse to infinite density lines 
and points can occur, collapse to an infinite density sheet is mathematically impossible. 
Correspondingly, in a two dimensional system, it is impossible to find collapse to an 
infinite density line, but it is still possible to observe collapse to an infinite density 
point [Hj. Let us briefly comment on this interpretation. Blow up in one dimension 
would mean that a real three dimensional system would be able to contract itself by 
using only one spatial dimension while the other two would remain invariant during this 
process. Global boundedness in time would invalidate this mechanism of chemotactic 
aggregation. Obviously, this fact crucially affects the patterns that can form. It is, 
however, important to point out that this connection between critical dimension and 
biological pattern formation is phenomenological, and a rigorous mathematical theory 
describing the interplay between dimensionality and geometry of the singular set for the 
Keller-Segel model constitutes an open problem [15]. So, assuming the phenomenology, 
the one dimensional result is not just a model problem, but has relevance in predicting 
the behaviour of the biological system. Interestingly, the three-dimensional modalities of 





dtc = DcAc + 7p — f3c. 
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chemotactic collapse allowed by the Keller-Segel system have already been observed in 
experiments performed with Escherichia coli [T7] . If this behaviour were general in 
cellular systems (collapse to points and lines, but not to sheets, in three dimensions, and 
to points, but not to lines, in two dimensions), and not a particular theory explaining 
these concrete patterns of Escherichia coli, we could name it as the Keller-Segel law 
for chemotactic aggregation. This would constitute an achievement of fundamental 
importance in mathematical biology. 

However, exceptions to this rule have already been found. In a recent in vitro 
experiment performed with mesenchymal cells [18] , it was shown that a two-dimensional 
system of these cells chemotactically aggregated into one- dimensional structures, clearly 
violating the Keller-Segel law. A theoretical explanation was found by taking into 
account the fact that diffusion is a too unrealistic assumption for mesenchymal cells [19] . 
Actually, these cells perform a sort of nonlocal diffusion that is better modelled by the 
operator A/(l — A), defined from its Fourier tranform 



instead of the Laplacian. Performing this substitution, it was proven that the 
corresponding nonlocal Keller-Segel system blows up in one dimension for large enough 
initial condition, leading thus to a modified Keller-Segel law that allows the type of 
chemotactic aggregation observed in the experiments [19] . 

At this point there is a very important question to be answered: is the Keller-Segel 
law a general result or a particular property of Escherichia coli colonies? It is clear 
that this law is valid whenever the Keller-Segel system applies, but what we would 
like to know is whether the Keller-Segel law extends beyond the validity of the partial 
differential system Eqs. (lla|16p . As happened in the case of mesenchymal cells, the 
diffusive approximation is too simplifying in many situations, particularly in those not 
designed at the laboratory. 

Diffusion is obtained as the description of the spatiotemporal distribution of a 
population density of random walkers [20]. However, in many situations found in 
Nature, Levy flights are commonly adopted as an efficient search strategy of living 
organisms [211 ESI EH]. Indeed, experimental evidence of superdiffusive behaviour has 
been found in biological systems [211 ESI [23]: super diffusion is characterized by a 
superlinear dependence in time of the mean square displacement of the position of the 
dispersing population. The correct description of a population undergoing Levy flights 
is given by the substitution of the Laplacian operator by a fractional operator of the 
Riesz type A"" in one dimension, defined from its Fourier transform 



where 1 < a < 2 ^24j. Direct substitution provides us in one dimension with the 
fractional Keller-Segel model: 




(2) 



(A"/)' 



(3) 





dtc = DcCxx + 7P - /5c. 



{Ah) 



The fractional Keller-Segel model 



4 



This model is thus presented as a more accurate description of a cellular population self- 
interacting chemotactically in those cases in which Brownian motion does not represent 
a good approximation for the population dispersal. So we can now translate our question 
to the present setting and ask ourselves if the Keller-Segel law is still valid for this partial 
differential system. It is a known mathematical fact that the Riesz operator A", with 
1 < a < 2 is less regularizing than the Laplacian, so it might be possible that the 
system Eqs.fl4a[ blows up in finite time. This would imply a new breakdown of the 
Keller-Segel law for an important number of situations. 

The goal of this article is to prove that solutions to a simplification of the system 
Eqs.fl4al l46l) are globally bounded in time. This is an important result since it implies 
the robustness of the Keller-Segel law for a large number of situations in which the 
original Keller-Segel system does not apply. The outline of the rest of the article is as 
follows. In Sec. |2] we present the simplified model and prove the simplest estimates. In 
Sec. [3] we prove the impossibility of finite time formation of Dirac masses, and in Sec. H] 
we prove that the solution is globally bounded in time. A summary of our results and 
some directions for future research are presented in Sec. [51 



tion. Changing 


variables 


X 




-(^) 
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iDp' 





2. Nondimensionalization and basic estimates 

The first step in the analysis of system Eqs. ( l4a[ |4M is to perform a nondimensionaliza- 

l/a 

X, (5a) 

{5b) 

p, (5c) 
and supressing the hats we arrive at the nondimensional fractional Keller-Segel model: 

dtp = AV - (pc^)x, (6a) 

^Jj^dtc = c..+p-^c. m 

Supposing that Dc » Dp (in the sense of much greater) and that n is large enough 
(but not making further assumptions on the value of (3), we finally arrive at the system: 

dtp = A"p - (pc^)x, (7a) 
Cxx = Sc- p, {7b) 



where 



This reduction to a parabolic-elliptic system is very common and has been performed 
many times in the literature [H O [25l [26], and it can be related to direct biological 
measures 
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To clarify our language, let us rigorously define what we understand by chemotactic 
collapse: 



Definition. We say that the system ( fTat 7b ) undergoes chemotactic collapse if 
there is a time < oo such that 

lim p(x,t) = +00, (9) 
for some x G M. 

The goal of this article is to prove that this system does not undergo chemotactic 
collapse. It is important at this point to note that there are other possible ways of 
defining chemotactic collapse [27], however, our definition is maybe more in the line of 

Refs. nain]. 

We take the domain to be the full line M, since boundary value problems for Levy 
flights are quite involved, as the long jumps make the definition of a boundary actually 
quite intricate (see however [28l[29]). Suppose that the initial condition po fulfils 

IIpoIIli < 00, (10) 

and as boundary conditions we take that p and c, and their corresponding spatial 
derivatives, vanish when |a;| 00. Since we are dealing with densities, we are only 
interested in nonnegative solutions to system Eqs. ( ITallTM . 

This system only models motion of the cells, but not any birth or death processes. 
This implies that the total number of cells, given by the norm of p, does not vary in 
time. It is actually easy to prove that the norm of p is conserved in time: 



Lemma. Let p, c be a solution to the system (7a, 7b). Suppose that for the initial 
condition ||po||li < OO; then 

IIpIIliW = IIpIIli(0), (lla) 

and 

l|c||Li(t) = I|c||li(o). (116) 
Proof. A direct application of the boundary conditions yields 

^llpllii = J A'^pdx- J (pc,),dx = 0, (12) 

because J A'^pdx = 0. Integrating Eq. (l76l) we find that 



|c||li = ^IIpIIli- (13) 



□ 



3. Impossibility of (5-function aggregates 



The aim of this section is to prove the following 

Theorem. Let p, c be a solution to the system 
condition obeys the properties stated in the last section and also fulfils ||p(a;, 0)11^2 < cxd. 



7a|, 7b). Suppose that the initial 
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Then there is a constant C depending on the and L} norms of the initial condition 
such that 

max{||p(-,t)|U2,||c(-,t)|U2}<C (14) 

for all times t > 0. 

Proof. Let us start evaluating the norm of p: 

Integration by parts yields 



ppxCxdx = I ppxCxdx + I p Cxxdx, (16) 
and using Eq.lff^ we see that 



— / pPxCxdx = [ p^cdx — - I p'^dx. (17) 



2 J ' 2 

We can use this result to find 

^|llp|li. = ilHli3-^ / 



-I|A"/Vlli.<^llp|li3-||A"/Vlli., (18) 

due to the positivity of c. The third moment of p might be estimated as follows 

llp|li3<llp|li^lHlLi, (19) 

now choose x ^ (1? Q^)? use a Sobolev embedding to find 

\\p\\u<cmh + \\^''^'p\\h)- (20) 



We can now invoke the Fourier transform of p 

1 



p{k) = I e'^''p{x)dx (21) 



to claim that 

l|A^/Vlli. = ll(A'^/V)1li. = / mp?dk = 
I \k\^p?dk+ I ^^\p\^dk< 

J\k\<R J\k\>R\k\''-'' ~ 

R''\\p\\l^ + ^ J \knp\'dk = R^\\p\\l. + ^\\A-/'p\\l., (22) 

where we have used the isometry of the Fourier transform in L^. We still need to 
estimate the second moment of p: 

Mlh < ||pIIli||pIIl- < ^IIpIIli + 

|llp|li^ < ^llpllii + c'- (iHii. + ||A"/VlliO ' (23) 
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where we have used a Sobolev embedding. Selecting e small enough we are led to 
conclude 



1 / 1 



2e' 

This inequality, in addition to Eg. (122!) yields 

\\p\\Is<c\\p\\l^ 



(24) 



;^IHIii + ^Jl|A"/Vlli.Vi5;^l|A"/Vlli. 



2e 



Ro'-x' 



(25) 



and substituting this result in Eq. (fT8|) we obtain 



d_ 

dt 



\p\\h < cWpWl^ 



(^Ilp|li.+^|||A"/Vlli.) + ^I|A"/V 



a/2^| |2 
L2 



2||A-/V||i2. 



(26) 



Now, by choosing a sufficiently large R and a sufficiently small e, we arrive at the 
estimate 



|iipiii.<c-ciiA"/viii., 

for suitable constants C, C > 0. Repeating the steps in Eq. fl23l) one finds 



-c|||A"'Vlli.<^l|p|li.-(i-c|)|H 



2 

L2' 



for small enough ei. Combining the last two equations yields 
d 



^llp|li.<c-c"||p| 



2 

L2; 



for suitable constants C, C" > 0, what in turn implies the desired estimate 

\\p\W<c. 

The boundedness of the norm of c can now be inferred from Eq. fl7&|) : 

5c = p+ c^x, 

what implies 



(27) 

(28) 

(29) 
(30) 
(31) 



|2 

Il2 



pcc/x + / CCxxdx < 



IIpI|l2||c||l2 - ||C3;.||^2 < ||p||L2||c||i:,2, 

and we finally have 

'IpIIl2 



C L2 < 



(32) 



(33) 



□ 



Corollary. This result precludes the formation of Dirac masses in the system 



m[76]). 
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We will prove in this section the global boundedness in time of the norm of p and c: 
Theorem. Let p, c be a solution to the system ( [7a| , 7b). Suppose that the initial 
condition obeys the properties stated in the last two sections and also fulfils 
dp{x, 0) 



dx 



< oo. 



(34) 



L2 



Then there is a constant C depending on the norm of the first spatial derivative of 
the initial condition, and the and norms of the initial condition such that 



max{p(a;, t), c(a;, t)} < C 

X 

for all times t > 0. 

Proof. Let us start estimating the time evolution of the norm of px 
1 d 



(35) 



2dt 



\px 



|2 

Il2 



PxdtPxdx = / pxA^pxdx 



/p.(pcJ..d. = -/p.(pcJ..d.-||A°'V.ll!.. 
The last integral in this equation may be rearranged in the following way 



(36) 



PxiP^x) xxdx 



PxPxx^^xdx 



2 / plcxxdx - / PxPCxxxdx = Ii + I2 + h 



(37) 

and we can manipulate each of the terms by reiteratively using Leibniz's rule, integration 
by parts, and the boundary conditions: 



h 
h 



plcdx 



25 / plcdx + 2 



plpdx, 



plpdx, 



and 



2 



p^cdx 



We finally find 



I = h + l2 + h 
6 



p^dx + / ppldx. 



(38) 
(39) 

(40) 



^2 

plpdx + / p'^cdx 



35 
Y 



p^cdx — - / p dx < - / p,j.pdx + 



52 



p^cdx, 



due to the positivity of p and c. Substituting in Eq. flS^ we get 



1 d 
2dt^ 



IP.IIi^ < I J Plpdx + ^ I p'cdx - \\A-/'px\\l.. 



It follows immediately that the second integral can be simply estimated by 
p'^cdx < ||p||i2||c||Loo, 



(41) 



(42) 



(43) 
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while the first integral might be estimated as follows 

plpdx < ||p^||L^||p^.||i2||p||i2 < ^\\p^\\l2 
We can use the Sobolev embedding 



h+2e\\p,\\% 



(44) 
(45) 



to arrive at 



- 2 

5 r 1 

2 2^ 



L2 



C LOO + 



(46) 



so Eq. fl36l) finally reads 



Ic?.. ^iiii9iiii ^ 



2 

L2 



+2eC(||p,.||i. + ||A"/V.||iO]-||A"/V.|li- 

We still need to estimate the norm of c. To this end we will use the 
transformed version of Eq. (l76|l : 

— k'^c = 6c — p, 

which implies 

This allows us to estimate the norm of c^ 



(47) 
Fourier 

(48) 
(49) 



|cx||l2 



k^\c\'^dk 



(P + 5)^ 



'dk < 



dk < c. 



(F + 5)2 

This fact, in addition to the Sobolev embedding 



lL°o 



< c 



Il2 



+ Ikajllia) , 



yields the desired inequality 

||c||l°o < c. 

This result, together with Eq.(l471). yields, for e small enough and suitable 
constants Ci, C2, and C3, the inequality 

|||p.||i.<Ci + C2||p.||i2-C3||A"/V.lli2. 
To finish our proof we still need the estimate 



(50) 

(51) 

(52) 
positive 

(53) 



\Mh= I \k\'p'dk= I \k\'p'dk + 



\k\<R 



\k\ 



2+a 



\k\>R 



\k\^ 



-p^dk < 



L2 + 



:I|A"/V.|li2, 



(54) 
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or, what is the same, 

-\\A'^/'p.\\h<R^'''\\p\\h-R''\Ml.. (55) 

By choosing R large enough, and combining this last inequality together with Eq.( l53il 
we obtain 

j^Wp.Wh < C, - C2\\p.\\l2, (56) 
for suitable positive constants Ci and C2, what implies 

\Ml^<C. (57) 
We recall again the Sobolev embedding 

llplli^ <c(i|p|ii2 + i|p.|iiO, (58) 

to conclude 

||p||loo<C. (59) 

□ 

Corollary. This result prohibits the chemotactic collapse for the system Eqs. 



m 7b) 



5. Conclusions 

In this work we have introduced and analyzed the fractional Keller-Segel model. This 
model is thought to be an extension of the standard Keller-Segel model to such situations 
in which the movements of the cells cannot be described by random walks, and a Levy 
flight description fits much better the cellular trajectories. Our analysis has revealed that 
solutions to the fractional Keller-Segel model exist globally in time in one dimension, as 
it happens with the standard Keller-Segel model, implying that chemotactic collapse to 
a sheet cannot occur in a three-dimensional setting, and chemotactic collapse to a line 
is impossible in two dimensions. This fact shows the robustness of the Keller-Segel law, 
that extends beyond its initial range of validity and applies to a much broader class of 
situations. 

We have performed our analysis taking the domain of the fractional Keller-Segel 
model to be the whole line M, which is a good first approximation (in the standard Keller- 
Segel model the proof of existence do not vary much from the unbounded domain to the 
bounded one f3]). However, it would be interesting to analize our model in a bounded 
domain, specially with Neumann boundary conditions, due to its biological relevance in 
experiments [3], [191 [26] . As shown in the diffusive case, homogeneous solutions become 
unstable for certain range of the parameter values, and a rich bifurcation structure, 
which may be supercritical or subcritical, arises [3]. 

Finally, we would like to remark that all the equations that have been written in 
this work are a continuous description of a discrete process. The number of cells involved 
must be integer, a fact that becomes crucially important when this number is not high 
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enough. This type of discrete processes admits a continuous description by means of 
stochastic partial differential equations [30], that favours its analytical assessment. It 
would be interesting to study the effect of such stochastic corrections on the problem of 
chemotactic collapse, to see whether or not they are able to modify the critical dimension 
of the problem. These and other considerations will be the object of future research. 
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